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Abstract. 

We study the structure of the phase diagram for systems consisting of 2- and 3- 
level particles dipolarly interacting with a 1-mode electromagnetic field, inside a cavity, 
paying particular attention to the case of a finite number of particles, and showing that 
the divergences that appear in other treatments are a consequence of the mathematical 
approximations employed and can be avoided by studying the system in an exact 
manner quantum-mechanically or via a catastrophe formalism with variational trial 
states that satisfy the symmetries of the appropriate Hamiltonians. 

These variational states give an excellent approximation not only to the exact 
quantum phase space, but also to the energy spectrum and the expectation values of 
the atomic and field operators. Furthermore, they allow for analytic expressions in 
many of the cases studied. We find the loci of the transitions in phase space from one 
phase to the other, and the order of the quantum phase transitions are determined 
explicitly for each of the configurations, with and without detuning. 

We also derive the critical exponents for the various systems, and the phase 
structure at the triple point present in the S-configuration of 3-level systems is studied. 


PACS numbers: 64.70.Tg, 42.50.Ct, 42.50.Nn, 03.65.Fd 
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1. Introduction 


While some observed phenomena such as the Rabi cycles in 2-state quantum systems 
may be explained by a semi-classical theory, other occurrences such as the revival of 
the atomic population inversion after its collapse [T]-[^ are quantum effects derived as 
a consequence of the discreteness of the field states. The revival property appears as 
well, for instance, in the dynamics of electron currents in monolayer graphene subject 
to a magnetic held [^. Even fractional revivals have been identihed with information 
entropies in different physical systems of interest [^. (A review of the formalism required 
to understand some aspects of the revival behaviour is presented in [^.) These and other 
purely quantum effects need to be studied through a quantum optics model such as the 
Jaynes-Cummings model (JCM) [^, which describes the behaviour of a 2-level system 
in the presence of a quantised radiation held. This model works very well when the 
radiation held and the system energy gap are close to one another and of the order 
of optical frequencies (~ 10^® Hz); this approximation is the so-called rotating wave 
approximation (RWA). The extension of the model to many “atoms” or systems is 
the Tavis-Cummings model (TCM) [^, and the removal of the RWA approximation 
including the so-called counter-rotating terms leads to the Dicke model (DM) , which 
describes the interaction of a single mode quantized radiation held with a sample of 
Na two-level atoms located inside an optical cavity, in the dipolar approximation (i.e., 
located within a distance smaller than the wavelength of the radiation). (Hereafter we 
will refer to “atoms”, but the theory applies to any hnite-level system, including spin 
systems and Bose-Einstein condensates.) The Dicke Hamiltonian has the expression 


H = hup a^a + uaJz ^ — {a^ J- + aJ+) -I —— {a^ Jp + aJJ). 


7 


( 1 ) 


Here, Na is the number of particles; the hrst term in the rhs represents the held 
Hamiltonian, where up is the held frequency and a\ a are the creation and annihilation 
photon operators; the second term represents the atomic Hamiltonian, with ua the 
atomic energy-level diherence, and the atomic relative population operator. The 2 
last terms represent the interaction Hamiltonian; we have written them separately in 
order to diherentiate the rotating term (hrst), with J± the atomic transition operators, 
and the counter-rotating term (second). 7 is the dipolar coupling constant. 

The parameters appearing in ([^ are related to the physical properties of the 
atom/system, and have dimensions as shown in Table where d is the dipole moment 
of the atom, e the electron charge, and p the atomic density inside the quantisation 
volume. It is convenient to redehne uja = 7 = and take up = 1 (i.e., measure 

frequency in units of the held frequency), which we do hereafter. 

The expression ([^ was derived from a multipolar expansion of the dipole interaction 
with the electromagnetic held. A diherent derivation related to the radiation gauge, 
where the long wavelength approximation is considered as well as the approximation 
to 2-level systems, leads to an extra diamagnetic term quadratic in the electromagnetic 
vector potential A, of the form -\-aY, with the diamagnetic coupling constant. 
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Table 1. Hamiltonian parameters in terms of physical parameters of the system. 


dimensional 


dimensionless 


'y = i^Ad\l^^ [freq] 'y = UAd ^ 


hujp 

2np 


hojp 


TT e^ p 


= o - KT = 

2m OJpJSa nbujpiy 

e^h r , e^h 


K = 


2m d?Cj\N^ 




[ 1 /freq] k = 


2m dP'CjAOJF Na 


This has led to some confusion in the literature as to the correct expression to use. 
Both the multipolar and the radiation Hamiltonians are related by a unitary gauge 
transformation, thus yielding the same physics; it is the approximation to 2 -level systems 
that breaks this symmetry (cf. for details). When using the Hamiltonian derived 
from the radiation gauge, the Thomas-Reiche-Kuhn sum rule would place contradictory 
bounds to the parameters of the model 11 ; furthermore, since the coupling strength 7 
is much smaller than the atomic level separation Oa for optical systems, it was believed 
that gauge invariance requires the presence of the diamagnetic term 12 . To the benefit 


of the Hamiltonian in Q, not only has a very strong case been made in its favour 
as a consistent description of the interaction of a one-mode light held with the internal 
excitation of atoms inside a cavity 13 , but experimental results indicate that transitions 
apparently forbidden by the no-go theorem from the sum-rule mentioned above can 
actually be observed 14,15 by using Raman transitions between ground states in an 
atomic ensemble. 

An important feature of atom-held interactions is the presence of phase 
transitions 16 from the normal to a collective behaviour: ehect involving all Na atoms 
in the sample, where the decay rate is proportional to N\ instead of Na (the expected 
result for independent atom emission). Quantum huctuations may drive a change in the 
ground state of a system, even at zero temperature, T = 0. A simple way to see this 
is to consider a Hamiltonian H{x), whose degrees of freedom vary as a function of a 
dimensionless coupling parameter y. The ground-state energy of H (y) would generally 
be a smooth, analytic function of y [^. Exceptions occur, for example, in the case 
when y couples only to a conserved quantity 


H{x) = Ho + xHi ( 2 ) 

where Hq and Hi commute. Then Hq and Hi can be simultaneously diagonalised, but 
while the eigenvalues vary with y, the eigenfunctions are independent of y. We can 
then have a level-crossing where an excited level becomes the ground state at a certain 
critical value of the coupling y = yc (cf. Fig. [^. 
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E 



E igure 1. Energy-level crossing 
state. 




For X ^ Xc and X ^ Xc the ground state of the system is clear; for x Xc terms 
start to compete, and the system would undergo a phase transition: we say that the 
limiting states realise distinct quantum phases. The crossing of levels in the spectrum of 
a physical system is an indication of a hrst order transition while, in general, the second 
order ones correspond to other causes (e.g., avoided crossings) and they are continuous. 
Each phase is a region of analyticity of the free energy per particle, and different phases 
are separated by separatrices which are singular loci of the free energy. Thus, the study 
of the phase diagram of a system is an important means to understand its behaviour. 

There have been various contributions to the study of phase transitions in 2-level 
In particular, the Husimi function has been used for phase space 


18-21 


systems 

analysis (22) and entropic uncertainty relations to detect quantum phase transitions 
Here, we want to stress the role of the catastrophe formalism to determine signihcant 
changes in the ground state of the system under small changes in the parameters of 
the model. Quantum phase transitions and stability properties have been extensively 
studied through the catastrophe formalism and the coherent states theory 


24 


31 


In 


particular, as these quantum systems cannot be solved analytically in an exact manner 
(except in the thermodynamic limit), in the latter references a procedure based on the 
use of the hdelity susceptibility of neighbouring states was established to determine with 


hne precision the location of the separatrices (but see also 19,3^). 

For applications such as quantum memories and other quantum information and 
quantum optics purposes it seems more appropriate to use 3-level atoms. Furthermore, 
approximations to 3-level systems in the A conhguration are plentiful: e.g., alkali metals, 
as conhrmed by the electromagnetically induced transparency effect. For practical 
applications, off-resonant systems protect one from spontaneous emission and have thus 
been favored because of their advantage when subjected to coherent manipulations; in 
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fact, schemes have been presented for various quantum gates using 3-level atoms and 
trapped ions 35,36 . The study of 3-level systems thus deserves attention. In particular, 


the importance of their phase diagrams has drawn the attention of some authors 33,34 


In 33 the energy surface method was applied to obtain an estimation of the ground 


state energy and the phase diagrams as well as the order of the phase transitions, in 
the three configurations, using the multipolar Hamiltonian. The results were compared 


with those of the exact quantum solution. In 34 the radiation Hamiltonian containing 


the diamagnetic term is used, in the Holstein-Primakoff realisation. They analyse the 
phase diagram of the three configurations in the thermodynamic limit, taking care of 
the regions where the Thomas-Reiche-Kuhn (TRK) sum rule holds, and they show that 
transitions from the normal to the collective regimes are possible even when the TRK 
rule is satisfied; this is in direct contrast to the situation of 2-level systems. 

Here, for the aforementioned reasons, we will consider the multipolar Hamiltonian 
and we will make use of the catastrophe formalism to study 3-level systems. Their 


Hamiltonian may be written as 37 

H = Hd + Hint, (3) 

where Hq and Hint are the diagonal and interaction contributions, respectively given 
by 

H£) = Vt of a -\- uji All + ^2 A 22 + 1^3 ^ 33 , (4) 


Hint 


—/== [fvi2 (^12 + A21) + fii3 (Hi 3 + ^31) + R23 {A23 + H32)] (a^ + a) 
vHa ^ ^ 


( 5 ) 


Here a\ a are as before the creation and annihilation electromagnetic field operators, 
and Aij = Mf fhe collective matter operators obeying the U (3) algebra 

[-djj, H-im] hji Aim him ^Iji (6) 

with a possible realisation Aif = and the total number of atoms is given by 


Na = X! 

k=l 


( 7 ) 


The Tth level frequency is denoted by Ui with the convention ui < U 2 < 003 , and the 
coupling parameter between levels i and j is fiij. The different atomic configurations 
are chosen by taking the appropriate value = 0 (cf. Fig. |^. 

We have written H (instead of ujp) for the frequency of the radiation field. The way 
in which equations (4]|5) are written lends itself to be easily generalised for a system of 
n-level atoms interacting with m-modes of a radiation field 

n n m 

-77general ~ ~ ^kj) (o'g + (8) 

i=l k=l j<k l=l V 

where the values of j, k, i are determined by the possible transitions according to the 
specific atomic configuration, and where we have Na = Yfk=i ^kk- 

In this work we shall review and extend the study of the phase diagrams for 2- and 
3-level quantum systems consisting of a finite number of atoms interacting through a 
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Figure 2. Atomic configurations S, A and V. The i-th level frequency is denoted by 
uji with the convention uji < UJ 2 < 1 ^ 3 , and the coupling parameter between levels i 
and j is gij. 


1-mode electromagnetic field. We show how the use of variational trial states that are 
adapted to the symmetry of the system Hamiltonian give an excellent approximation 
not only to the exact quantum phase space, but also to the energy spectrum and the 
expectation values of the atomic and held operators. When in the RWA approximation, 
the total number of excitations is an integral of motion of the system; using trial states 
adapted to the symmetry of the Hamiltonian then means essentially projecting onto this 
integral of motion. In the full model (rotating and counter-rotating terms), however, it 
is the parity in the number of excitations that is conserved, and to obtain symmetry- 
adapted states (SAS) we therefore take linear combinations of coherent states of the 
same parity. 

These symmetry adapted states were hrst used in 


38 , named “even and odd 


coherent states”, as nonclassical states for the study of singular non-stationary quantum- 
mechanical harmonic oscillators, and later to discuss the properties of the tomographic 
representation of quantum mechanics 39,40 . Here, we use them to look in detail at the 


structure of the phase diagram and the behaviour of the phase changes. We also present 
some virtues and limitations of these symmetry-adapted states, use the hdelity and the 
hdelity susceptibility of neighbouring quantum states to hnd the loci of the transitions 
in phase space from one phase to the other, and derive the critical exponents for the 
various systems. 


This work is dedicated with great appreciation to Professors Vladimir and Margarita 
Man’ko in their joint 150-year celebration, for their numerous contributions to the 
development and promotion of guantum optics and mathematical physics. 


2. Two-level Systems 

The simplest completely soluble quantum-mechanical model of one 2-level atom in an 
electromagnetic held is described by the Jaynes-Cummings (JCM) model (^. This, and 
its generalisation to Na identical 2-level atoms given by Tavis and Cummings |^, the 
TCM model, were fundamental to study basic properties of quantum electrodynamics 
and to understand phenomena like the existence of collapse and revivals in the Rabi 
oscillations (observed experimentally for the hrst time in 1987 [^). Both the JCM and 
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the TCM models discard the terms in the Hamiltonian which do not conserve the total 
number of excitations of the held plus matter by using the RWA approximation. When 
these terms are considered we obtain the full Dicke model (DM) [^. In this Section we 
consider the phase diagrams presented by these models for the ground state, both in 
the case of a hnite number of atoms and in the thermodynamic limit, by making 
use of the catastrophe formalism to determine when signihcant changes to the ground 
state occur for small changes of the external environment (the parameters of the model). 
The inhuence of the phase transitions on the behaviour of observables of interest for the 
matter and the held are also presented. 

The choice of the use of the catastrophe formalism allows us to obtain analytic 
descriptions for the phase diagram in parameter space, which distinguishes the normal 
and collective regions, and which gives us all the quantum phase transitions of the ground 
state from one region to the other as we vary the interaction parameters (the matter-held 
coupling constants) of the model, in functional form. This approach thus allows also for 
the study of the asymptotic behaviour in any of the quantities of interest: the number 
of particles, the constants of motion, and the interaction parameters themselves. 

Catastrophe theory derives from the research of Rene Thom in topology and 
diherential analysis on the structural stability of diherentiable maps 41 . Dissipative 


systems, for example, always reach equilibrium; this equilibrium is characterised by a 
certain function fi{x) which at x represents the minimum of usually the energy of the 
system, and when this minimum fi{x) is stable x will be a regular point in the space of 
parameters describing the system. But when the energy changes abruptly at fi{x) due 
to slight disturbances the local minimum is destroyed in a neighbourhood of x, p,{x) 
ceases to be an attractor of the dynamics, and x is a catastrophic point: the state of the 
system will present sudden jumps from x to another point x' (another attractor) and 
back. The dynamics of the system thus bifurcates. It is these bifurcations that we are 
interested in studying analytically. 


2.1. The Jaynes-Cummings and the Tavis-Cummings Models 

A 2-level system of atoms interacting dipolarly with an electromagnetic held of 
frequency up is described by the Tavis-Cummings Hamiltonian [^, which we may write 
as 

where we have set h = 1 and all quantities are dimensionless. We have also divided the 
expression by Na in order to consider an intrinsic Hamiltonian, which we do hereafter. 
We can set cup = 1 (i.e., measure frequency in units of the held frequency), and dehne 
a detuning parameter A = up — ua = 7 — oj a] thus, A = 0 when particles and held 
are in resonance and A 7 ^ 0 when away from resonance. It is convenient to introduce 
A = \J.P -\- 1/4 — 1/2 -|- J 2 -|- a^a because it turns out to be an integral of motion for 
the system. Its eigenvalues are A = u m j, with j = A/2, j + m the number of 
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Figure 3. Colour online. Energy spectrum for the Hamiltonian Q of the Tavis- 
Cummings model. A ^ 0 (left) leads to avoided crossings. Pairs of curves of the same 
colour emanating from almost the same point on the energy axis correspond to the 
same value of A. When A = 0 (right) curves corresponding to the same A touch at 
7 = 0. The straight horizontal line at E = —0.4 (magenta) is the energy of the ground 
state in the normal region. 


atoms in their excited state, and v the number of photons. The Hamiltonian can then 
be rewritten as 


H = 


-A 


A 


Jz + 


7 


(a^ J_ + aJjf). 


( 10 ) 


Na" Na‘ ' s/N2Na 
The eigenvectors and eigenvalues of H can be obtained through diagonalisation 
of its associated matrix, thus allowing us to calculate the expectation value of all 
important field and matter observables, as well as the entanglement entropy, the 
squeezing parameter, and the population distributions 25,26 . For instance, taking the 


natural Hilbert space basis |z/, j, m), where u is the eigenvalue associated to the photon 
number operator, j{j + 1 ) is the eigenvalue associated to the total angular momentum 
operator, m is the particle occupation number \m\ < j < Na/‘2 , where the j = Na/2 
holds for identical atoms. Substituting the label m for the eigenvalue of the constant 
of motion, A = z/ + m + j, we can obtain the full energy spectrum of H. This is shown 
in Figure (left) for Na = 6 atoms, A up to 10, and a detuning parameter of A = 0 . 2 . 
One can see the avoided crossings due to A 7 ^ 0 ; had we A = 0 they would touch at 
7 = 0 (cf. Figure 1^ (right)). Pairs of curves of the same colour emanating from almost 
the same point on the energy axis (or the same one in the case for A = 0) correspond 
to the same value of A. The thicker horizontal line at = —0.4 (magenta) is the energy 
of the ground state in the normal region. 

We are interested, however, in studying the system analytically. To this end, we 
propose to use as a test-state a direct product of coherent Heisenberg-Weyl HW{1) 
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states and SU{2) states |a,C) = |a) ® |C) as: 

I ^ exp(-|a|V 2 ) “ ^ ^ / 2 j 

(1 + ^0 Jp'-j W' + 


1/2 

^j+m 


v) ® |j, m) 


( 11 ) 


with C = tan exp(i0), a = + ip), (6^,0) being the parameters on the Bloch 

sphere and {q,p) the held qnadratnres. 

The energy surface, dehned as the expectation valne of the Hamiltonian on the test 
state: PL = {a, C\H\a, (), is then given by 


11 'y 

+1^^) “ 2^^ COS0 + sin 0 (geos 0-psin0). (12) 

The critical points of Pi determine 3 regions, as given by 0c = 0 {North Pole), 
9c = 71 {South Pole), and 9c = arccos(a;A/7^) {Parallels)] for each of these regions the 
minima of the energy Eo and valnes Ac := (A)c (the expectation valnes of the constant 
of motion) are as follows: 


9c = 

0, 

Eo = - 

Na 

2 ’ 

Ac = 

0, 

9c = 

TT, 

Eq = 

Na^a 

2 ’ 

Ac = 

Na, 

9c = 

arccos 

, Eq — - 

NaVI+P) 

472 ’ 

Ac = 

77a(—<^A ('^A+2)-I-7^-I-27^) 

472 


for uja > 7^ 
for ua < -7^ (13) 
for Icn^l < 7^ 


At these critical points, Qc = — \lN a/'^ 7 sin 9c cos 0c, and pc = ^ Na/2 7 sin 9c sin 0c, 
so that matter and held variables combine. As 0 is a cyclic variable, 0c may be taken 
arbitrarily. We set 0c = 0, and the expressions in terms of this variable may be recov¬ 
ered by performing a rotation throngh an angle 0 aronnd the z-axis in the appropriate 
phase space: {q, p) and {J^, Jy) for held and matter qnantities respectively. 

We can write explicitly the form that the states take in each of these 3 regions: 

North Pole: ua > 7^ iV’np) = |0) ® \j, -j) 

Sonth Pole: ua < - 7 ^ If^sp) = |0) ® \j, j) (14) 

Parallels: I^aI < 7^ If’par) = Emi-/ Y.t=o^rn,u W) ® \j, m) 


with 


A 


m,i/ — 



1/2 

exp{ 


J7^ 

4 



(j-m+iy)/2 



(-v/2j7) 

Vfd 

(j+m+iy)/2 


The 3 regions dehne also a separatrix, where the Hessian of PL is singnlar. This 
is given by uja = 4=7^, and is shown in Fig. Crossing the separatrix along paths 
I, II, III, and IV (horizontal (green) and vertical (brown) straight lines in the hgnre) 
leads to second-order phase transitions; crossing it along path V to hrst order transitions. 

In general, these coherent variational states approximate very well the properties 
of the gronnd state of the qnantnm solntion 25 . This is trne for the energy, the 


constant of motion A( 7 ), and the matter observables {Jf) and its flnctnation sqnared 
(AJ^)^, etc. Even the expectation valne of the nnmber of photons n = (Nph) is well 
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Separatrix 

0)A 



Figure 4. Colour online. Phase diagram for the Tavis-Cummings model. The normal 
region is described by the North and South Poles, separated by parabolae from the 
Parallels which denote the collective region. These parabolae constitute the separatrix. 
Crossings along paths J, II, III, and IV (horizontal (green) and vertical (brown) 
straight lines) lead to second-order phase transitions; only the crossing V (slanted 
(orange) straight line) through the origin gives a first-order phase transition. 


approximated; but its fluctuation, as well as other properties of the system such as the 
occupation probabilities, are not: Fig. (left) shows how bad an approximation to the 
photon number fluctuation we get. The noticeable differences arise from the fact that 
the coherent state contains contributions from all eigenvalues A = i/-hm-hjofA, and 
therefore does not reflect the symmetry of the Hamiltonian leading to the constant of 
motion. 

One may maintain the symmetry through a projection of the variational tensorial 
product of coherent states onto the value of the constant of motion of the TCM which 
minimises the (classical) energy of the ground state. This projection restores the 
Hamiltonian symmetry and is amiable to analytical calculations. 

Projecting, the state becomes 


|0) 0 jj, -j) , 


= J\f< 

I ^v= 


,] (a-,)' \p A - J - I/) , IcuaI < 7^ 


uja > 7 


i/=max[0, X-2j] \\- 


(15) 


|0) ® \j, j) , wa < -7^ 

where we have defined p = — F^'>' (i -|- ^). xhe factor Af is the normalisation factor. 
With respect to these projected states, the energy surface is given in terms of associated 
Laguerre polynomials 26 as follows 


n 

H 


A - j + j A _ 

2j 



27 

ysj’'. 


. i<A<2i 

, A > 2j 


(16) 


A = 0 


(17) 
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{^nf|NA 



(An)VA'A 



Figure 5. Colour online. Left: Comparison between the fluctuation in the expectation 
value of the number of photons (An)^/ACi given by the coherent state approximation 
(continuous, red curve) and the exact quantum solution (discontinuous, blue curve), 
shown as functions of 7 . Right: The same comparison using projected states 
(continuous, red line); note the scale in the ordinate axis. We have used Na = 20, 
and A = 0.2 in both graphs. The noticeable differences arise from the fact that the 
coherent state contains contributions from all eigenvalues A = + m + j of A, and 

therefore does not reflect the symmetry of the Hamiltonian. 


and the approximation to the photon nnmber flnctnation is restored as shown in Fig. 
(right) (note the scale of the ordinate axis). 

A better way to measnre the “distance” between states is via the hdelity, 


F (pi, Q 2 ) = tr Q 2 \/^ 


where qi and Q 2 denote the density matrices of the states in qnestion. For pnre states, 


this dehnition coincides with the sqnare of the scalar prodnct between the states 19 


Fignre [271] shows a perfect overlap F = 1 between the projected and qnantnm states in 
the normal region, dropping to F = 0.996 when crossing the separatrix into the Parallels 
region, only to recover again towards F = 1 as 7 grows. 

Even if onr approximation by projected (symmetry-adapted) states is not exact, 
an excellent approximation to the exact qnantnm solntion of the ground state of the 
TCM model is obtained. What is gained is that these states have an analytical form in 
terms of the model parameters and allow for the analytical calculation of the expectation 
values of field and matter observables, as well as for the study of the phase diagram of 
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Figure 6. Fidelity F between the projected state and exact quantum ground states, as 
a function of the interaction strength 7 . The plot is for a detuning parameter A = 0.2 
and Na = 20 atoms. 


the system. 


2.2. The Dicke Model 


When the RWA approximation is not taken, we have the full Dicke Hamiltonian given 
in equation ([^. Once again, one may obtain analytical expressions for the energy and 
expectation values of the relevant operators of the system via the use of the Heisenberg- 
Weyl and SU{2) coherent states (0 as trial states, and the variational procedure 
described above. This trial state contains N = 2j particles distributed in all the possible 
ways between the two levels and up to an inhnite number of photons in the cavity. The 
energy surface in this case is given by 

and the separatrix shrinks to ua = for 7 ^ the critical value of 7 . As before, 
the crossings of this separatrix are second-order phase transitions, except for the hrst- 
order crossing through the origin. The energy minima in the normal and collective 
(superradiant) regions are 


Enormal 2 ,Na'^i 

,2 


E superradiant Ahl7 

and the expected number of photons are 

{-h^ph) normal 0 

{■^ph) superradiant A^TT 


F - 


- T 


(19) 


( 20 ) 


which calls for the dehnition of a: = y/yc. 











Phase Diagrams of Systems of 2 and 3 levels in the presence of a Radiation Field 13 


(An)VA'A r 




Figure 7. Colour online. Left; comparison between the fluctuation in the expectation 
value of the number of photons (An)^/iVyi given by the symmetry-adapted state 
approximation (continuous, red curve) and the exact quantum solution (discontinuous, 
blue curve), shown as functions of 7 . We have used Na = 20, and A = 0. Right: F 
functional dependance on x for different values of Na- As Na increases F tends to 
zero much more rapidly. 


AsA = yj2 + l/4-l/2+J,+ata is no longer a constant of motion for the system we 
cannot simply project onto one of its eigenvalues, rather, we have a dynamical symmetry 
associated with the projectors of the symmetric and antisymmetric representations of 
the cyclic group C 2 , given by 
1 


= 


(1 ± e*'"') 


( 21 ) 


This symmetry allows, however, for the classihcation of the eigenstates in terms of the 
parity of the eigenvalues \ = j + m + voi K 


27 


Adapting the coherent states to the 
parity symmetry of the Hamiltonian then amounts to sum over A even or odd, with two 
resulting orthogonal states |q;, C, ±). For these states the energy surface associated to 
the superradiant regime takes the form 

-4n 


{H)± = -N^y 




1±D 


with 


jr — 


'(l-x '‘j 


( 22 ) 


(23) 


and the limit a; —)■ 1 gives the expressions for the normal region. The hdelity between 
these symmetry-adapted states and the exact quantum states is very close to 1 except 
in a small vicinity of the transition region in phase space, so it is no surprise that they 
provide an excellent agreement with the expectation values of the quantum operators 
for the system, an example of which is shown in Fig. (left) for the fluctuation in 
the expectation value of the number of photons (Anj^/Ajq as given by our projected 
state approximation (continuous, red curve) compared with the exact quantum solution 
(discontinuous, blue curve), as functions of 7 . We have used Na = 20, and the resonant 
condition A = 0. 

If we calculate the overlap between the coherent and adapted states we obtain 

KocC I ocC, 3 =) r ^ 2 d 1=.^) 


(24) 
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and, since the behavionr of D falls very rapidly with 7 (cf. Fig. (right)), this 
overlap will be at best eqnal to 1 / 2 , which makes the ordinary coherent states a good 
approximation only in special cases. [Appendix A compares the expectation valnes and 
flnctnations of matter and field observables for the coherent and symmetry-adapted 
states, evalnated at the critical points for the energy snrface (18). For expectation 


valnes different from zero in the symmetry-adapted states, the coherent state results 
can be obtained from the former by letting D go to zero. Notable exceptions are the 
field quadratures (g, p) and the atomic operator fluctuations. For large Na the function 
T tends to zero even more rapidly; this is why coherent states have been so successful 
in the past as trial functions. 

Like the quantum states, the symmetry-adapted states show no divergences for field 
or matter expectation values at the phase transition. This is in contrast with results 
found previously 15,21,42 , which are an artifact of an inappropriate truncation of 


the Hamiltonian. For more good properties of the symmetry-adapted states, including 
probability distributions of photons, of excited atoms, and their joint distribution, 
cf. (^. In particular, even though the coherent states, the symmetry-adapted states, 
and the quantum states, are quantities arrived at via very different methods, they show 
a universal character in that a universal parametric curve for any number of atoms Na 
is obtained for the first quadrature of the electromagnetic field, q, and for the atomic 
relative population (Jz), as implicit functions of the atom-field coupling parameter 7 , 
valid for both the ground- and first-excited states 29 . Furthermore, for all values of 


the coupling parameter and again any number of atoms, the behaviour of the number 
of photons vs. the relative atomic population is universal. 


2.2.1. Critical Exponents For a homogeneous function /(r) we have f{(3r) = g{(3) /(r) 
for all values of (3. The scaling function g{/3) is of the form g{(3) = (3^; s is called the 
critical exponent. It is known that the singular part of many potentials in physics are 
homogeneous functions near second-order phase transitions; in particular, this is true 
for all thermodynamic potentials 43 . The behaviour of important observables of a 


system near phase transitions may thus be described by the system’s critical exponent, 
and these are believed to be universal with respect to physical systems. 

Our treatment for finite 2-level systems in a cavity, in the presence of a radiation 
field, allows us to study the critical value of the atom-field coupling parameter 7 ^ 
as a function of the number of atoms Na, from which its critical exponent may be 
derived. Figure shows this relationship for the ground state of both the quantum 
states (left) and the symmetry-adapted states (SAS) (right). For the quantum states 
the points correspond to a numerical solution from diagonalising the Hamiltonian, and 
the continuous curve to a model fit. The value of 7 c was obtained by calculating in 
parameter space the place where the fidelity function between neighbouring states (cf. 
equation ([30)) below) vanishes. We varied Na from 10 to 1800 and the logarithm of the 
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Figure 8. Logarithmic behaviour of the critical value of the coupling parameter 7 c 
with the number of atoms Na- For the quantum (q) states the critical exponent is 
—2/3, while for the symmetry-adapted states (SAS) it turns out to be —11/21. 


variables is plotted for a more demanding fit, obtaining 


In 


Tc 


= In 


In(iVA) , 


or, eqnivalently. 


7/ = - + A iv 
2 2 ' 


1 -2 
' "a ^ 


(26) 


(26) 


Except for a small vicinity of the phase transition, the SAS states do approximate 
very well the qnantnm solntions. However, the critical exponent obtained for the 
asymptotic behavionr of the adapted states is —11/21, as opposed to —2/3, as shown 
in the hgnre (right). This is precisely becanse the evalnation takes place at the phase 
transition point, where the states (qnantnm and adapted) differ most 10 . The valne of 


7c for the SAS states was obtained by calculating in parameter space the place where the 
minimum of the energy E+min for the state |acCc, +) presents a discontinuity. Since we 
are interested in the asymptotic behaviour, we took Na from 200 to 1000; the continuous 
curve shows the £t 


7c 


= 2+"' 


■5 Na 


(27) 


Table 1^ shows a sample of values of {Na, 7 c) for the quantum and the SAS ground 
states, in order to make explicit the fact that for small Na the values of the quantum 
critical interaction parameter 7/ differ considerably from that of the SAS states 7/“®. 
This difference tends to zero as Na increases, and in the limit Na —>■ C)0 the phase 
transition region in phase space coincides for both states at 7c = 0.5. 

It is interesting to note, from equation (24), that only for 7 = 0 (i.e., no matter-field 
interaction) do we have 


(oc q I «£ t. ±) t = 1 (7 = 0). 


(28) 


i.e., the overlap between coherent and symmetry-adapted states is perfect only when 
the interaction Hamiltonian Hint vanishes. As soon as there is an interaction, no matter 
how small, the states differ. This is due to the fact that the ground states coincide only 
at 7 = 0. Even in the normal regime, where the coherent ground state has exactly zero 
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Table 2. Sample of values of the quantum critical interaction parameter 7 ® and the 
SAS critical interaction parameter 7 ®“^, for different values of Na- 


Na 

7/ 

^sds 

1 c 

20 

0.5677 

0.5522 

40 

0.5432 

0.5343 

100 

0.5236 

0.5204 

400 

0.5096 

0.5097 

800 

0.5061 

0.5068 

1000 

0.5051 

0.5060 



Figure 9. Energy per particle (left) and the expected number of photons per particle 
(right) for the symmetry-adapted ground state inside the normal region. The (red, 
dashed) straight line at EAmin/NA = —0.5 is a comparative reference for the energy 
of the ground coherent state. We have taken Na = 10. 


photons, the SAS ground state (just as the quantum ground state) is a superposition of 
states with an expectation value for Np^ different from zero. This is true for any hnite 
number of atoms Na- Figure]^ shows the energy per particle and the expected number 
of photons per particle for the ground symmetry-adapted state inside the normal region. 
We have taken Na = 10 to make the distinction visually clear. 

In the asymptotic limit a: —)■ 00, equation (24) gives a value of 1/2 for the overlap of 
the coherent and the SAS ground states. The same is true in the limit Na —>■ C)0. This 
is to be expected, as the SAS ground state has contributions only from the even-parity 
components of the coherent ground state. 


3. Three-level Systems 

A 3-level system of Na atoms interacting dipolarly with an electromagnetic held of 
frequency D is described by the intrinsic Hamiltonian given in equations (|^|^|^. Once 
again, we may take 0 = 1 and measure all frequencies in units of the held frequency. As 
mentioned before, the Ath level atomic frequency is denoted by uji with the convention 
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1^1 <1^2 <1^3, and the coupling parameter between levels i and j is The three 
different atomic conhgurations are chosen by taking the appropriate value /ijj = 0 (cf. 
Fig.g). It is also convenient to define a detuning parameter = Ui — Uj — D between 
levels i and j. 


(29) 


In the RWA approximation the Hamiltonian reduces to 37 

H = Vt a^a + uji An + uj2 A22 + 1^3 A33 
1 

and it has 2 constants of motion, viz., the total number of atoms Na = ^u-, and 

the total number of excitations M = afa + A2A22 + A3A33, where the value of A, {i = 2, 3) 
depends on the configuration taken (cf. Table [^. 


/ii2 (a A21 + A12) + 7^13 (n A31 + ^13^ + (123 (a A32 + A23) 


Table 3. Values of A^, z = 2, 3, for the constant of motion M in the different 
configurations. 


Configuration 

A2 

-^3 

- 

1 

2 

A 

0 

1 

V 

1 

1 


Notice that the Hamiltonian (29) is invariant under the transformation a —)■ —a 
and af —?• —a\ which preserves the commutation relations of the bosonic operators. For 
this reason we consider only positive values for /ijj. As the system cannot be solved 
analytically, one may solve via numerical diagonalization. A natural basis in which we 
diagonalize our Hamiltonian is \v] q, r) 33 . Here, v represents the number of photons 


of the Fock state; r, q — r and Na — q are the atomic population of levels 1, 2, 3, 
respectively. 

In order to study the phase diagram of the system we make use of the fidelity F 
and the fidelity susceptibility y of neighboring states 32,44 , defined by 

F(r, t + 6t) = |('0(r)|'0(r + 5r))|^ 


_ o 1 ^ 

^ (5r)2 

Whereas the hdelity is a measure of the distance between states which vary as functions 
of a control parameter r, the fidelity susceptibility, essentially its second derivative with 
respect to the control parameter, is a more sensitive quantity. The fidelity measure goes 
to zero at each phase transition, as the nature of the ground state changes completely 
and orthogonaly; the fidelity susceptibility has divergences at these critical points in 
phase space. Crossing a separatrix produces a change in the total excitation number 
(M). 
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To follow a similar procedure as for the 2-level systems, and be able to study the 
phase diagram analytically, we consider as a variational trial state the direct product 
of Heisenberg-Weyl HW{1) coherent states for the radiation part, |a} = |0), and 

U (3) coherent states constructed by taking the exponential of the lowering generators 
acting on the highest weight states of U{2>) 

1C} := I[^i,^2,/i3]7n72,73} = [^i,^2,/i3]) , (31) 

where | [hi, ^ 2 , hs]) represents the highest weight state of the Gelfand-Tsetlin basis 
in the irreducible representation [hi,h 2 ,h 3 ] of U{3) |^. When we consider the 
totally symmetric representation [N^, 0 , 0 ] (for indistinguishable particles) the trial state 
becomes 

|a; C) = |«; W,72,73} = e““'| 0 ) ® |[W, 0 , 0 ]) (32) 

where the parameter 71 no longer appears since ^32 | [Na, 0 , 0 ]) = 0 . 

It is convenient to use a polar form for the complex parameters 


a-.= p exp(i 0) , 7 j := pj exp(i (pj) , j = 1, 2, 3 , (33) 

and minimising with respect to these new parameters the energy surface 'H{a, () = 
{a; Cl H \a; C}/{«; C|ctj C} fhe RWA approximation takes the form 


1 


^rwa(Pc) P2ci PSc) — T { 

i\a 


iVi 022 p\c + ^^3 pic 




P 


P12 P2c + Pl3 P3c + P23 P2c P3c }/ (l + pL + pL) ) 


(34) 


where pc, P 2 c and p^c denote the critical values of the corresponding variables, and we 
have taken pi = 1. It is important to stress that equation (34) is valid for all three 
conhgurations. 

From this minimal surface the hrst separatrix corresponding to the phase change 


M = 0 —)■ M 7 ^ 0 (i.e., from the normal to the collective regimes) is given by 33 

i) for the S-con£guration 

P12 + [ IP23I ~ ^/^31 ] 0 (IP23I ~ = ^21 ; 

ii) for the A-con£guration 

P13 + [ IP23I ~ \/^2l ] 0 (IP23I ~ \/^2l) = a ^31 ; 

iii) for the R-configuration 

P12 P13 _ ^ . 

U21 UJ31 

where ujij = Ui — ujj and 0 is the Heaviside function. These are shown in Fig. 10 
for double atomic resonance with respect to the radiation held frequency in the S- 
conhguration (i.e., 0231 = 2 0221 = 2 ), a small detuning in the A-conhguration (0221 = 
0.2, a; 3 i = 1), and double atomic resonance in the R-con£guration (0221 = 0231 = 1). For 
equal atomic detuning in the A-con£guration the separatrix is identical to that of the 
R-con£gur at ion. 


(35) 

(36) 

(37) 
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A - config 

^23 





^0-2, A^12 - A^23 “ 0.2 

X 



Figure 10. Shape of the separatrix M = 0—J-MT^Oin phase space for the different 
configurations. We have taken double atomic resonance, except for the A-configuration 
where we took a detuning of 0.2 between one of the transitions and the field frequency. 
The triple point in the S-configuration is shown (see text). The fidelity susceptibility 
for neighbouring states (lower right) is shown for the (projected) SAS states in the 
S-configuration along the path /212 = ^23 — 0.2. 


Further analysis shows that in the S-conhguration the phase transition across 
hi 2 = yhRfi is of second-order, while that across the segment of the circumference 
is of hrst-order. In the A-con£guration with unequal atomic detuning we have the 
same behaviour. In the A-configuration with equal atomic detuning and in the Id- 
configuration, however, all transitions are of second-order. 
























Phase Diagrams of Systems of 2 and 3 levels in the presence of a Radiation Field 


20 



Figure 11. Fidelity between the quantum and (projected) SAS ground states (left), 
and that between the quantum and coherent ground states (right), for Na = 3, V- 
configuration. 


Being (M) = v + \ 2 {^ 22 ) + ^ 3 (^ 433 ) a constant of motion, we obtain states adapted 
to the symmetry of the Hamiltonian by projecting onto the appropriate value of (M). 
This is done in practice by substituting z/ = (M) — \ 2 {A 22 ) — ^ 3 (^ 33 ) and keeping the only 
relevant value of (M). These are the (projected) SAS states in the RWA approximation. 

In the thermodynamic limit, given by z/ oc Na with Na —)■ 00 , the loci in parameter 
space of the quantum phase transitions are exactly those shown by Fig. But even for a 


small number of atoms the approximation to the separatrices given by the projected SAS 
states is remarkably good: the hgure shows, in its lower right, the hdelity susceptibility 
divergences at each phase crossing along the path /ri 2 = /r 23 — 0.2 for the S-con£guration 
and Na = 2 , as a function of p.23- Since /ii 2 = 1 is hxed and independent of Na at this 
separatrix {vide infra), the projected SAS prediction gives /i23 = 1-2 which compares 
well with the value of (123 = 1-28 for the hrst transition of the exact quantum ground 
state, even though Na = 2 . This good approximation by the chosen variational states 
obeys the fact that the hdelity between the quantum and projected SAS ground states 
gives a perfect overlap except in a small vicinity of the phase transitions, as shown in 
Figure [IT| (left). The reader may compare this vs. the overlap between the quantum and 
the coherent ground states shown at right. In both cases Na = 3 and we have chosen 
to illustrate the result for the R-conhguration (those for the other conhgurations being 
very similar). 

All phase transitions tend to those given by equations (35, 36, 37) as Na —)■ 00 
with z/ oc Na- In this thermodynamic limit these are the only ones that remain. For the 
I/-con£guration the consequent transitions take place at a family of curves congruent 
with and ever more distant to the one shown in Fig. IT, which approach the latter 
uniformly as Na grows. For the S-conhguration we have curves with similar shape to 
that shown in Fig. 10, with a vertical straight edge and an upper circular arc. The 
vertical edge tends to that at /ii2 = 1 as Na grows, while the circular arcs “slide down” 
the /i 23 -axis, intersect the arc of the transition M = 0 —)■ M 7 ^ 0 shown, and continues 
sliding down this arc tending to (/ii2, 1^23) = (1, a/2) as Na —)■ 00 . We show this in 
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Figure 12. Critical values gijqc of the interaction parameters for the quantum 
transitions M —)■ M + 1 as functions of the number of atoms in the S-configuration. 
This shows how the loci of the quantum phase transitions change as the number of 
atoms grow. In the limit Na —> oo they converge to the separatrix between the normal 
and collective regions. 


Figure 12 The subfigure at left shows the critical value /ii 2 qc of /ri 2 for the quantum 


transition M —)■ M + 1 as a function of the number of atoms, i.e., how the transitions 
to the right of the straight vertical line /ii 2 = 1 move as Na changes. They all tend to 


the limit fii 2 qc = 1, as given by equation (35) when Na —)■ oo. At right we plot p 23 


i qc 


as a function of Na, to see how the phase transitions above the circular arc move; the 
first transition isM = 0—)-M = 2 since the phase region M = 1 stops at /i 23 = 
and does not reach the upper arc. We see that the point where these phase regions 
intersect the circular arc slide towards /i 23 = \/2, again as given by equation (35). In 


the thermodynamic limit, then, the separatrix reduces to the line segment given by 
fii 2 = 1 and /i 23 € [0, ^/2], plus the arc of circumference starting at /i 23 = \/2. 

The A-configuration has a similar behaviour as the Id-configuration when in double 
resonance, and a behaviour much like that of the S-configuration when away from double 


resonance. 


3.1. A Triple Point in Phase Space 

The S-configuration is special in that it shows a richer structure. In particular, it has a 
triple point in parameter space, corresponding to the place where the phases for M = 0, 
M = 1, and M = 2 meet 47 . The term triple point is mainly used in the context of 
fluids, where different phases of the fluid meet in parameter space. Here we use the same 
terminology since the different values of the total excitation number M correspond to 
completely different structures of the ground state, even though the energy is the same 
for all of them, and since these three regions meet at a point in parameter space as shown 
in Figure 13 The meaning is also the same as for a thermodynamic triple point: any 


fluctuation (in this case quantum) will drastically change the composition of the ground 
state. And since in the collective region we have a decay rate proportional to Na^, as 
opposed to Na for the normal region, this gives hope for experimental exploitation of 
the triple point. 
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In the RWA approximation and in double resonance this triple point resides at 


(/^i 2 , ^ 23 ) = (1, \/2) (ch Fig. 10). In the full model, contemplating the counter-rotating 


terms, we just divide these values by 2 {vide infra, Subsection 3.2). It is a fixed point, 
independent of Na, which subsists in the thermodynamic limit. It is also characteristic 
of the S-configuration; it does not appear in the A or the V configurations. We can 
calculate the ground state |' 0 )c/s at the triple point for each phase, by diagonalising 
the Hamiltonian in the basis |z/; q, r). In the RWA approximation one gets analytic 
expressions. For Na > 2 we have: 

M=0: 


I gs 


= | 0 ; Na, Na) 


(38 a) 


M=l: 


gs 


M=2: 


I gs 




N^. Na) 


(386) 


| 1 ; Na, Na 


1 ) + 


1 Na-1 


Na 


|0; Na, Na — 2) -|- 


t) + 2 1^’ 


Na) 


(38c) 


It is clear that, even when we are at the same point in phase space, the ground state may 
acquire very different structures. Away from double resonance the triple point is still 
present, though its coordinates in phase space vary as well as the specific combination 
given by the equations above. 

When the number of excitations M is small the dimension of the Hilbert space does 
not depend on Na, making it possible to study the system in the limit Na —t 00 . The 
energy spectrum, in particular, does not depend on (123 in this limit, and it shows a 


collapse of energy levels at precisely /ii 2 = 1 for all values of M. Figure shows this 
for M = 0 to M = 5, and it is interesting to compare it with the spectrum of the 2-level 
Tavis-Cummings model. Fig. As a function of M, at the triple point, we have an 
equidistant spectrum with only even harmonics 47 , and it is interesting to note that 


at /ri 2 = 1 we have precisely all the even harmonics as degenerate energy levels, and no 
others. 

The behaviour at the axis /ii 2 = 0 is also interesting. The total degeneracy for each 
M found at /ri 2 = 0 in the limit Na —t 00 only survives for finite Na at p .23 = 0, i.e., 
when there is absolutely uo matter-field coupling. As soon as the coupling is “turned 


on”, this degeneracy breaks down. Figure 15 shows E vs. /X 12 for Na = 4 (left) and 
Na = 100 (right), when /i 23 = 0 (blue, coutiuuous line) and when 1 x 23 = 1-5 (red, dashed 
line). While the ground and first excited states still show degeneracies at /ii 2 = 0, these 
are broken for the second excited state. 
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a M=o 

□ M=1 
■ M=2 


Figure 13. Energy of the ground state plotted as a function of y .12 and y ,23 for the 
5-configuration in double resonance. The 3 regions meet at a point, the triple point, 
at coordinates (l,\/2, 0 ) in parameter space (marked in the figure with a black dot). 
This point is independent of the number of atoms, and subsists in the thermodynamic 
limit. 


E 

10 

8 

6 

4 
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Figure 14. Colour online. Energy spectrum in the limit —>■ 00 of the 5- 

configuration. A collapse of energy levels for all values of M at precisely ^12 = 1, the 
surviving separatrix in phase space, is clear. Different colours correspond to different 
values of M; equally, the value of M can be read as the value for the energy at y,i 2 = 0. 
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Figure 15. Colour online. Energy spectrum degeneracies at y,i 2 = 0 broken when 
Na is finite. Shown are Na = 4 (left) and Na = 100 (right), when ^23 = 0 (blue, 
continuous line) and when /i 23 = 1.5 (red, dashed line). 


3.2. Counter-Rotating Terms: the full model 


When we do not make the rotating wave approximation, i.e., we include the 
counter-rotating terms in equations Hill, minimising the energy surface PL^a, Q = 
{a; Cl H \a] C}/{«; C|ctj C} with respect to the polar parameters takes the form 


'^(Pc, P2c, P3c) = + { 


UJi + UJ2 pIc + ^3 pIc 




pc 


P12 P2c + Pl3 P3c + P-23 P2cP3c }/ (l + pL + pL) ) 


(39) 


Comparing equations (39, 34), the energy surfaces Pi and PL^na coincide if we 
identify 


(P'i^)RWA ^ ^ (Pjfc) ■ (40) 

This means that "Hrwa will inherit the properties of Pi at values of (pij )bwa equal to 
\pij. (This is the same behaviour as that mentioned earlier for the Dicke model.) In 
particular, the shape of the phase diagram will be inherited in full at coordinates half 
those of the RWA scenario, and the order of the phase transitions will be the same. 

Whereas M is a constant of motion in the RWA approximation, it is not in the full 
model. As in the 2-level DM model, it is the parity in the number (M) of excitations 
that is conserved, as U{9) := exp (iPM) is only a symmetry operator for P = 0, tt. 
To obtain symmetry-adapted states we therefore take linear combinations of coherent 
states of the same parity 


l«; C}± := (1 ±exp[i7rM]) |a; C} , 


( 41 ) 


and the energy surface for these SAS states, in any conhguration, results in 46 




± 


= — D |a|^ [exp(|aH (7* • 7)^^ =F exp(-|aH (7* • 7)^-^ 

to! A 

+ 2 ^Ui |7i|^ exp(|aH (7* • 7)^^^"^ ± (-1)^' exp(-|a|^) (7* • 7)^^"^ 

i=l 
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i<j=i 


X 


\Na — ^ 


(42) 


exp(|a| ) (7* + 7* 7i) (7* • 7) 

± exp(-|«n ((-!)"» 7; 7^. + (-1)"^' 7; 7,) (7* . 7)^^-i' 

where 7 = (71, 72, 73), 7 = (7i, (- 1 )'^'' 72 , (-1)'^®73), and 71 = 1. Again, one may 
nse the polar form of these parameters to minimise with respect to each one in order 
to obtain the minimnm energy snrface for the system. In general this has to be done 
nnmerically, bnt the 17-con£gnration lends itself to an analytic treatment; fnrthermore, 
all transitions in this conhgnration are of second order, making it a good candidate for 
the stndy of its critical exponents. 


3.3. Critical Exponents 


Using the polar form given in eqnation (33), and fnrther dehning 


p 2 =i cos(r 7 ), = i sin( 7 ), 

P 12 = P cos(6'), /ii 3 = p sin(6'), 
X =r]-0 


(43) 


the SAS energy snrface (42) in the U-con£gnration takes the form 


Wv + - 


(i-eCh+e) i-e+Pp-i+e)) 

(-1 + e) (1 + eC («' +ppi+e)-ipip cos(x) 

where we have dehned 

,Na 


1C = 




+ e 


2NaP‘ 


(i+U 


Na 


(44) 


(45) 


We minimise with respect to the parameters (p, f), and having an analytical expression 
for the SAS gronnd state allows to evalnate the relevant held and atomic operators in 
this state. Of particnlar interest is the comparative behavionr of the system in the 
normal and collective regimes: Fignre [l6| (left) shows the (normalised) atomic nnmber 
of excitations in comparison with the (normalised) field excitations. This snggests that, 
in the normal region, 

(A22 + ^ 33 )-!- _ {Nph)+ 

(An)+ Na 

as is to be expected from the atomic decay rate in a non-collective regime; as soon as the 
system enters a collective regime the nnmber of held excitations increases mnch more 
rapidly than the atomic excitations. Differently to the coherent states, the SAS and 
qnantnm gronnd states in the normal regime contain non-zero contribntions of photonic 


(46) 


and atomic excitations. Eqnation (46) thus may be used, in general, as a criterion 


to define the normal region by using the appropriate atomic operators for the atomic 
excitations in each conhgnration. 
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Figure 16. Colour online. Left: equality of the (normalised) atomic and field 
excitations in the normal regime; upper (blue) curve corresponds to the field excitations 
and lower (red) curve to the atomic excitations. Right: discontinuity in the phase 
parameters at a phase transition for the symmetry-adapted ground state in the V- 
configuration; upper (blue) curve corresponds to Pj^in lower (red) curve to Cmin- 


A phase transition canses a change in the strnctnre of the gronnd state, which is 
reflected by a discontinnity in the phase parameters (see Figure (right)). We use 
this discontinuity to find the critical value of the interaction strength /i^ at the phase 
transition that separates the normal from the collective regimes as a function of the 
number of atoms, from Na = 100 to Na = 2000. Figure [T^ shows a logarithmic plot for 
these two variables, together with the linear £t 




sas 

c 


1 

2 


3 —^ 

-b e 2 21 


(47) 


It is interesting to compare this relation with that obtained for the Dicke model 
using SAS states, equation (27). The critical exponent is exactly the same. An analysis 
of residuals shows a conhdence interval of [—0.530, —0.519] for the exponent, to a 
conhdence level of 0.95, where —11 fits perfectly, with a goodness-of-fit = 0.9997. 


4. Discussion and Conclusions 

We have reviewed and expanded upon the structure of the phase diagram for systems 
consisting of 2- and 3-level particles dipolarly interacting with a 1-mode electromagnetic 
field, inside a cavity, paying particular attention to the case of a finite number Na 
of particles, and showing that the divergences that appear in other treatments are a 
consequence of the mathematical approximations employed, and can be avoided by 
studying the system in an exact manner quantum-mechanically or via a catastrophe 
formalism with variational trial states that satisfy the symmetries of the appropriate 
Hamiltonians. 

We have shown how the use of these variational states give an excellent 
approximation not only to the exact quantum phase space, but also to the energy 
spectrum and the expectation values of the atomic and field operators. Furthermore, 
they allow for analytic expressions in many of the cases studied, even for finite Na- We 
have made use of the fidelity and the fidelity susceptibility of neighbouring quantum 
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Figure 17. Logarithmic behaviour of the critical value of the coupling parameter 
/ic with the number of atoms iV^, for the symmetry-adapted ground state in the V- 
configuration. The linear fit (continuous straight line) shows a critical exponent of 
— 11/21, exactly as that found for the Dicke model using SAS states. 


states to find the loci of the transitions in phase space from one phase to the other; 
having analytic expressions allows for the order of the qnantnm phase transitions to be 
determined explicitly for each of the conhgurations, with and without detuning. Finally, 
we have derived the critical exponents for the various systems. 

The S-configuration in 3-level systems is particular in that it exhibits a triple point 
in phase space. This means that any quantum fluctuation at this location will drastically 
change the composition of the ground state. The exact form of the ground state at this 
triple point has been studied; the same can be done for excited states in the vicinity of 
this point or elsewhere in parameter space. 

Finally, a criterion (equation (46)) to define the normal region for the full 
Hamiltonian in the different conhgurations was suggested, acknowledging the fact that 
the SAS and quantum ground states in the normal regime contain non-zero contributions 
of photonic and atomic excitations. 

With the promise of the benehts of quantum information the study of these systems 
acquire greater importance, as they constitute the basic q-dit blocks themselves as well 
as the possible quantum logical gates for computational purposes. The properties of the 
systems treated here have been intriguing, no less because of the search for a hne control 
of the light-matter interaction at the level of single and few atom-photon pairings. It is 
hoped that this manuscript conveys an accurate account of the properties and structure 
of the phase space of these interesting systems. 
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Appendix A. Expectation values and fluctuations: Dicke Model 


The following table shows the expectation values and fluctuations of matter and 
field observables for the coherent and symmetry-adapted states in the collective 
(superradiant) regime. The mean-field behaviour for the normal region can be recovered 
by taking the limit a; —)■ 1. (Adapted from (48|.) 


Comparing equation (12) and equation (18) it is clear that, by substituting 
7c -A- 7c/2 in the Table, we obtain the expectation values of operators for the Tavis- 
Cummings model, and their fluctuations. 




Table Al. Expectation values and fluctuations of matter and field observables for the coherent and symmetry-adapted states in the 
superradiant regime. The mean-field behaviour obtained in the normal region can be recovered by taking the limit a; —)■ 1. 
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